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ABSTRACT 

We present the results from the analysis of the broad-band X-ray spectra of 5 Anomalous 
X-ray Pulsars (AXPs) and Soft 7 -ray Repeaters (SGRs). We fit their Suzaku and INTEGRAL 
spectra with models appropriate for the X-ray emission from the accretion flow onto a pulsar. 
We And that their X-ray spectra can be well described with this model. In particular we And 
that: (a) the radius of the accretion column is ~ 150-350m resulting in a transverse optical 
depth of ~ 1; (b) the vertical Thompson optical depth is 50-400, and (c) their luminosity 
translates in accretion rates 10'^gs“'. These results are in good agreement with the predic¬ 
tions from the fail-back disk model, providing further support in the interpretation of AXPs 
and SGRs as accreting pulsars. 

Key words: pulsars: individual (4U 0142 h-61, IE 1547-5408, 4U 0142h- 61, IRXS J1708- 
4009, SGR 1900H-14) - X-rays: stars - stars: magnetic fields - accretion disks 


1 INTRODUCTION 

The anomalous X-ray pulsars (AXPs) are objects that show X-ray 
pulsations but, in comparison to the more luminous accreting pul¬ 
sars in X-ray binaries, they exhibit several unusual properties (e.g. 
no evidence for a companion star, softer X-ray spectra than ac¬ 
creting pulsars). Soft 7 -ray Repeaters (SGRs) on the other hand 
are sources that produce short (< 1 sec, typically ~ 0.1 s) bursts 
with luminosities up to ~ 10 "" erg s^'in low-energy 7 -rays and 
X-rays. They also exhibit rare, giant flares with luminosities up 
to 10 “*“* erg s“', as well as low-luminosity ( 10 ^^^^® erg s“*) persis¬ 
tent emission (e.g. Woods & Thompson 2006). It is now widely 
accepted that SGRs and AXPs are manifestations of the same class 
of objects (e.g. Kouveliotou et al. 1998; Woods & Thompson 2006) 
which, in addition to the bursting activity, share some general com¬ 
mon characteristics. These common traits include: low X-ray lu¬ 
minosities (10^"*^^® erg s“*), but significant hard X-ray emission 
(e.g. Meregheti 2008; Enoto et al. 2010), relatively long periods 
(P~ 1-12s), and large spin-down rates (P> 10^'^ss^*) (e.g. Woods 
& Thompson 2006). 

These properties have been generally interpreted within the 
context of the magnetar model (e.g. Duncan & Thompson 1992; 
Thompson & Duncan 1995; Beloborodov 2013). According to this 
model, AXPs are together with SGRs, two manifestations of non¬ 
accreting young pulsars with ultra-strong magnetic fields in the 
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10"^^’®G range (e.g. Woods & Thompson 2006, and references 
therein). The persistent emission of magnetars is believed to be due 
to the decay of their strong magnetic field, while their outbursts are 
interpreted as violent magnetic reconnection events due to the re¬ 
arrangement of plates on the crust of the pulsar. The ultra-strong 
magnetic field of these systems is inferred from their spin decay. 
However, recently several members of the AXP/SGR family have 
been found to have dipole magnetic fields (as derived from their P 
and P) in the more typical range of ~ lO’^^'^G for accreting pul¬ 
sars (e.g. Rea et al. 2014; Rea et al. 2013). This has led to the notion 
that the activity of these objects is powered by an internal toroidal 
field, which also leads to the formation of local multipole fields. 
Although the magnetar model explains the energy output and the 
outbursts of AXPs and SGRs, their broad-band X-ray spectra are 
usually described on the basis of ad-hoc phenomenological mod¬ 
els such as combination of power-law components (e.g. Meregheti 
2008; den Hartog et al. 2008a, 2008b; Enoto et al. 2010). How¬ 
ever, so far there has not been any physical model that reproduces 
consistently the entire broad-band X-ray spectra of AXPs/SGRs. 
For example, Beloborodov (2013) proposed a spectral model based 
on the outflow of e“-e^ pairs produced close to the surface of 
the neutron star, and it was recently applied to the spectra of 3 
AXPs/SGRs (Hascoet et al. 2014), providing a good description 
of their hard X-ray spectra, but requiring the ad-hoc inclusion of 
multiple modified black-body components to fit the spectrum be¬ 
low lOkeV. On the other hand the “twisted magnetosphere” model 
of Zane et al. (2009), while it fits well the spectrum below 10 keV, it 
requires an additional power-law component to reproduce the hard 
part of the spectrum. 
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An alternative model to that of an isolated spinning-down 
magnetar is that of a newborn pulsar with a magnetic dipole field 
(IqH-iSg)^ accreting leftover gas either from the common-envelope 
phase (van Paradijs et al. 1995), or from the fall-back disk of its par¬ 
ent supernova explosion (Chatterjee et al. 2000; Alpar 2001). In this 
case, the spin down of the neutron star is attributed to the interaction 
between the accretion disk and the magnetosphere (e.g. Ertan et 
al. 2012). Owing to the low luminosities, the accretion rate of these 
systems is expected to be much lower than that of typical accreting 
pulsars in X-ray binaries. In addition, given their low luminosity 
and lack of donor stars, we would not expect to observe strong op¬ 
tical emission from these systems, in agreement with the general 
characteristics of the AXP/SGR population. In this context, small 
scale outbursts observed in AXPs/SGRs can be easily explained 
in terms of accretion of clumps in the fail-back disk. On the other 
hand, large outbursts are explained in the context of reconnection of 
local multipole magnetic fields anchored on the surface of the neu¬ 
tron star. In general this “accretion model” reproduces well the spin 
period distribution, light curves, and broad-band multi-wavelength 
emission of AXPs and SGRs (e.g. Ertan et al. 2007; ^ah§kan & 
Ertan, 2012). 

The discovery of hard X-ray emission from AXPs/SGRs ex¬ 
tending up to > 100 keV (e.g. Kuiper et al. 2004, 2006; den Hartog 
et al. 2008; Enoto et al. 2010), gives us an additional test-bed for 
the fall-back disk model. In the first of a series of papers exploring 
the consistency of the hard X-ray emission from AXPs/SGRs with 
the “accretion model” (Trumper et al. 2010), we found that indeed 
the hard X-ray spectra of the brightest AXP, 4U 0142-1-61, can be 
interpreted in terms of accretion. The Chandra and INTEGRAL 
spectrum of this source can be well represented by a model includ¬ 
ing thermal as well as bulk-motion Comptonization (bmc model; 
e.g. Mastichiadis & Kylafis 1992; Titarchuk et al. 2007). A follow¬ 
up study of the energy-dependent pulse profiles of the same source 
(Trumper et al. 2013) showed that they can be well modeled in the 
context of the fan beam from the accretion column, and a polar 
beam from the hot polar cap. More recently, Guo et al. (2015) an¬ 
alyzed a sample of 4 AXPs/SGRs with the same spectral model as 
in Trumper et al. (2010). They found that the hard X-ray spectra 
of these 4 additional sources can also be well reproduced with the 
bmc model, further supporting the accretion model for AXPs. 

However, the spectral models used in the studies of Trumper 
et al. (2010) and Guo et al. (2015) have several shortcomings when 
applied to accreting pulsars: they assume spherical or disk geome¬ 
tries, instead of the cylindrical geometry of the accretion flow onto 
a pulsar, and they do not account for the different scattering cross- 
sections in the presence of a strong magnetic field. In order to fur¬ 
ther explore the consistency of the fail-back disk model with the X- 
ray spectra of AXPs and SGRs we embarked on a systematic study 
of their broad-band X-ray spectra, taking advantage of the current 
generation of hard X-ray telescopes {Suzaku and INTEGRAL ), and 
employing the best available spectral models for the X-ray spec¬ 
tra from the accretion flow onto an accreting pulsar. In this paper 
we present the spectral analysis of the 0.5-200 keV spectra of 5 
AXPs/SGRs with spectral models appropriate for modeling the ac¬ 
cretion flow onto a pulsar. 

The structure of this paper is as follows: in §2 we present the 
sample and data used in this study; in §3 we describe the analysis 
of these data, and the results from the analysis of the X-ray spectra; 
and in §4 we discuss the implications of our results for the accre¬ 
tion model of pulsars. Finally in §5 we summarize our findings. 
All quoted errors correspond to the 90% confidence limit for one 
interesting parameter. 


2 THE SAMPLE AND DATA 

In order to obtain useful constraints on the X-ray spectral prop¬ 
erties of AXPs/SGRs we selected to study objects from the list 
of known AXPs/SGRs with available X-ray spectra above 10.0 
keV (e.g. Mereghetti 2008; Enoto et al. 2010). We searched the 
Suzaku (Mitsuda et al. 2007) and INTEGRAL (Winkler et al. 2003) 
archives for available data. We found 5 objects with good quality 
spectra above lOkeV. These are listed in Table [T] along with the 
basic information of the observations used in this study. One addi¬ 
tional object with known hard X-ray emission (SGR1806-20) was 
excluded from our analysis because it is near the Galactic Center 
and its hard X-ray emission might be significantly contaminated by 
the Galactic Ridge emission (Krivonos et al. 2007). 


2.1 Suzaku data 

The Suzaku payload includes two X-ray detectors: the X-ray Imag¬ 
ing Spectrometer (XIS; Koyama et al. 2007), and the Hard X-ray 
Detector (HXD; Takahashi et al. 2007, Kokubun et al. 2007). The 
XIS consists of four identical CCD-based instruments (three of 
which are operational), which provide imaging and spectroscopic 
data in the 0.5-10 keV band. It is placed at the focal point of the X- 
ray Telescope, which provides focused images with a half-power 
diameter of 2' (Koyama et al. 2007). The HXD consists of two in¬ 
struments: a PIN silicon-diode detector sensitive in the 10-60 keV 
band, and a GSO/BGO phoswich scintilator sensitive above 40 keV. 
Being a non-imaging instrument, it uses collimators in order to re¬ 
duce the background. The lack of imaging information complicates 
the analysis of data for sources close to the Galactic Ridge since its 
contribution to the spectrum of the sources can only be modeled. 

The Suzaku data were obtained from the HEASARC archive. 
For their analysis we followed the standard procedures described 
in the Suzaku ABC guidJH v. 5.0, using the FTOOLS v.6.11 data 
analysis suite. The XIS data were first reprocessed and screened us¬ 
ing the aepipeline FTOOL with the default parameters. Then 
we extracted lightcurves for each of the 3 available XIS detectors 
(XISO, XIS2, XI3) in order to search for strong background flares. 
None of the sources showed any significant flaring that would re¬ 
quire additional screening. 

Then we extracted images for each of the 3 XIS detectors 
based on which we defined source and background apertures. The 
source apertures were selected to include as many as possible of the 
source photons, and in general their radius was ~ 3 arcminutes. Lo¬ 
cal background was selected from circular regions adjacent to the 
source, while taking care to avoid spillover of source counts due to 
the wings of the Point Spread Function (PSF). From those regions 
we extracted the source and background spectra used in our analy¬ 
sis. Response matrices for the source spectra were calculated using 
the xisrmfgen and xissimarfgen FTOOL. 

The HXD data were reprocessed to apply the latest calibration 
files, and screened using the aepipeline FTOOL. The standard 
screening criteria listed in the Suzaku ABC guide were applied. 
Since the HXD is not an imaging instrument, an estimate of the X- 
ray and particle backgrounds is based on detailed modeling of the 
in-orbit background data (Fukazawa et al. 2009). The final back¬ 
ground spectrum was calculated from the non-X-ray background 
data produced by ISAS for each observation sequence, and a model 
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Table 1. Sample and data used in this study 


Object 

RA Dec 

Suzaku 

INTEGRAL 



(J2000) 

Obs Date Exposure 

Obs. Period 

Exposure 


hh mm ss.s dd mm ss.s 

XIS, PIN 





(ksec) 


(Msec) 


IE 1547-5408 

15 50 54.1 

-54 18 23.8 

2009-01-28 

10.6, 33.4 



IE 1841-045 

18 41 19.3 

-04 56 11.2 

2006-04-19 

97.9, 63.4 

2003-03-2009-11 

2.341 

4U 0142+61 

01 46 22.4 

+61 45 03.3 

2009-08-12 

107.4, 99.7 

2002-12-2008-04 

1.35 

IRXS 11708-4009 

17 08 46.9 

-40 08 52.4 

2009-08-23 

60.9, 47.9 

2003-02-2010-03 

2.7 

SGR 1900+14 

19 07 14.3 

+09 19 20.1 

2006-04-01 

21.7, 14.3 

2003-03-2010-03 

3.1 


Positions from Mereghetti (2008); the exposure times are corrected for deadtime. 


of the cosmic X-ray background (Moretti et al. 2009) using the 
hxdpinxbpi FTOOL. 


2.2 INTEGRAL data 

The International Gamma-Ray Astrophysics Laboratory (INTE¬ 
GRAL) carries the Imager on Board the INTEGRAL Satellite 
(IBIS; Ubertini et al. 2003) and the Spectrometer on INTEGRAL 
(SPI; Vedrenne et al. 2003). These instruments use coded aperture 
masks allowing the reconstruction of X-ray images in the hard X- 
ray band. The IBIS instrument in particular consists of a CdTe de¬ 
tector array (ISGRI; Lebrun et al. 2003) and a Csl detector array 
(PICsIT). The ISGRI detector that will be used in this study, is sen¬ 
sitive in the 20 keV - 1 MeV energy range, and provides a field of 
view of 19° X 19°with a spatial resolution of ~ 12'. The ISGRI 
detector is more sensitive than PICsIT, but overall it is less sensi¬ 
tive than the Suzaku PIN detector that operates in their overlapping 
energy range. 

We obtained the INTEGRAL ISGRI spectra and response ma¬ 
trices from the HEAVENS database (Walter et al. 2010), which 
provides data products for known sources from the analysis of all 
publically available INTEGRAL data up to the time of download 
(in our case September 2013). However, since these data are taken 
over the course of several years and generally at a different epoch 
from the Suzaku data, there is always the possibility for spectral 
or intensity variations. Eor this reason we generally prefer to use 
the Suzaku XIS and HXD data. When there is agreement between 
the spectra of Suzaku XIS and INTEGRAL ISGRI we use both in 
our analysis in order to improve the statistics of our spectra and to 
benefit from the extended coverage to higher energies provided by 
ISGRI. 


3 SPECTRAL FITS 

3.1 Description of the physical model 

Since, according to the fail-back disk model, the energy output 
of AXPs/SGRs is generated by accretion of material onto the 
pulsar, one would expect X-ray spectra similar to those of X- 
ray pulsars. The expected accretion rates from a fail-back disk 
(~ 10*^ - lO'^g s“*; Ertan et al. 2009) are in the range of accretion 
rates of low-luminosity accreting pulsars. In this case, one would 
expect a power-law spectrum with a high-energy cutoff, which re¬ 
sults from thermal and bulk-motion Comptonization in the accre¬ 
tion flow (e.g. Lyubarskii & Sunyaev 1982; Becker & Wolff 2007, 
hereafter BW07; Kylafis et al. 2014). The slope of the power law 


and the energy of the cutoff depend on the energy of the seed pho¬ 
tons, the temperature of the electrons, the velocity of the accretion 
flow, and the optical depth in the accretion column, which deter¬ 
mines the number of scatterings that a photon will undergo before 
emerging (e.g. BW07; Kylafis et al. 2014). 

The spectral model required to fit the data corresponds to the 
following physical picture (first discussed in Triimper et al. 2010, 
2013; Kylafis et al. 2014): photons produced in the accretion 
flow are up-scattered by the thermal as well as the bulk-motion 
Comptonization mechanisms by the infalling electrons. These up- 
scattered photons comprise the power-law component of the ob¬ 
served spectrum. In addition to the non-thermal emission from the 
accretion flow, thermal photons from the polar cap of the neutron 
star also contribute to the observed spectrum, either directly, or af¬ 
ter beeing comptonized by electrons at the atmosphere of the polar 
cap (c.f. Eerrigno et al. 2009). 

Modeling the X-ray spectra produced by accretion of matter 
onto a pulsar is a very complex and still open problem. The best 
currently available models for the spectra produced in the accre¬ 
tion flow are based on the treatment of Becker & Wolff (Becker 
& Wolff 2005; BW07) which was implemented in XSPEC by Eer¬ 
rigno et al. (2009). This model used the Green’s function formalism 
to solve analytically the radiative transfer equation within a cylin¬ 
drical accretion flow in the presence of a strong magnetic field. 

As discussed in BW07, the photons produced at or below the 
accretion shock are trapped in the accretion flow and they can only 
escape sideways. The optical depth in the transverse direction is 
of the order of a few (e.g. BW07; Triimper et al. 2013; Kylafis 
et al. 2014), therefore the photons escaping from the sides of the 
accretion column have undergone several scatterings and they have 
increased their energy. This results in a hard power-law spectral 
component associated with the fan beam. On the other hand, the 
cross-section along the magnetic field axis is (BW07) 

/ GTiE/Ecf if E<Ec 
ar if E>Ec 

where Ec = 1L 6 S 12 keV is the cyclotron energy, and 612 is the mag¬ 
netic field at the surface of the neutron star in units of lO'^ G. Since 
the magnetic field along the magnetic axis decreases (as B oc z^^), 
the energy of the resonant cyclotron scattering cross-section will 
decrease accordingly. This makes the vertical direction of the ac¬ 
cretion flow (i.e. parallel to the dipole magnetic field axis) opti¬ 
cally thick to all photons produced at the thermal mound or below 
the accretion shock. 

The Eerrigno et al. (2009) model includes approximations for 
the angle-averaged scattering cross-sections in the parallel and per¬ 
pendicular directions with respect to the magnetic field, and it ac¬ 
counts also for the effects of thermal and bulk-motion Comptoniza- 
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tion (e.g. Titarchuk et al. 1997), the latter being particularly im¬ 
portant in this accretion flow geometry. It accounts for the three 
main emission mechanisms that operate within the accretion flow: 
thermal (black-body) emission at the base of the accretion column 
(thermal mound), cyclotron emission due to transitions from the 
first Landau level of the collisionally excited electrons in the strong 
magnetic field of the pulsar (also including the inverse process of 
resonant cyclotron absorption), and bremsstrahlung radiation from 
the electrons. However, it has several shortcomings such as decou¬ 
pled treatment of the thermal and bulk Comptonization, various ap¬ 
proximations for the relation between parameters of the radiative 
transfer and the geometry of the accretion flow (e.g. parameteriza¬ 
tion of the velocity of the accretion flow in terms of the optical 
depth), angle-averaged scattering magnetic cross sections, and a 
simple cylindrical accretion flow. Despite these shortcomings, until 
recently, this was the best available model for analyzing the spectra 
of accreting pulsars. 

More recently, Farinelli et al. (2012; hereafter F12), following 
the basic formulation of the radiative transfer equation by BW07, 
proposed another method to calculate the spectra from the accre¬ 
tion column, that overcomes some of the limitations of the Ferrigno 
et al. (2009) implementation. More particularly, this new solution 
allows for a more flexible treatment of the velocity profile of the ac¬ 
cretion flow by allowing the velocity to be a function of the height 
rather than a function of the optical depth, and most importantly 
reduces the number of fitted model parameters that are strongly 
correlated. 

Therefore, in our parameterization of the spectrum, we com¬ 
bine the F12 model which describes the emission from the accre¬ 
tion column, and a black-body component which describes the di¬ 
rect thermal emission from the polar cap. One would expect that a 
fraction of the photons emitted by the polar cap would interact with 
the accretion flow and be Comptonized, contributing to the hard X- 
ray emission. However, as discussed in Kylafis et al. (2014), for a 
cylindrical geometry of the accretion flow, the fraction of these iter- 
sected photons is very small and they hardly contribute in the hard 
X-ray spectra, supporting the decoupled modeling of the thermal 
emission from the polar cap and the emission from the accretion 
column that is adopted in this study. 

The F12 model is parameterized in terms of several physical 
parameters: the temperature of black-body photons at the thermal 
mound at the base of the accretion flow (kT™), the velocity profile 
of the accretion flow, the temperature of the electrons in the ac¬ 
cretion flow responsible for the thermal Comptonization {kTe), the 
total vertical optical depth along the accretion column (i.e. paral¬ 
lel to the magnetic field lines; r), and the radius of the accretion 
column (ro). 

For the velocity profile, the F12 model allows for two possi¬ 
bilities. One is the velocity profile adopted by BW07, which is 

^(r) = -'I'r, (1) 

where {3 is the velocity of the infalling material in terms of the 
speed of light, VE' is a proportionality factor, and r is the verti¬ 
cal optical depth. This means that the magnitude of the velocity 
is maximum at the top of the accretion column, i.e. at the place 
which BW07 call “sonic surface” (see their Fig. 1), and zero at the 
surface of the neutron star. Such a velocity profile is appropriate for 
high-luminosity X-ray pulsars, in which the radiative shock occurs 
far from the neutron star surface. In this model, the flow above the 
“sonic surface” is not taken into account. 

The second velocity profile allowed in the F12 model is the 
free-fall one: 


/3(z) = -A(-r, ( 2 ) 

z 

where z ^ Zs is the distance from the center of the neutron star, Zs is 
the distance of the sonic point in the accretion flow from the cen¬ 
ter of the neutron star (c.f. Basko & Sunyaev 1976), rj is the index 
of the velocity profile (rj = 0.5 for free fall), and A is a normaliz¬ 
ing constant defined as A = PoiR/zs^, where /3o is the magnitude 
that the velocity, in units of the speed of light, would aquire at the 
surface of the neutron star of radius R. Such a velocity profile is ap¬ 
propriate for low-luminosity X-ray pulsars, in which the radiative 
shock occurs near the neutron star surface. In this model, the flow 
below the sonic point is not taken into account. 

In our analysis, we use the second velocity profile, because it 
is more appropriate for the low luminosities involved in AXPs and 
SGRs. 

Another parameter of the F12 model is the albedo of the sur¬ 
face at the base of the accretion flow, which for a neutron star is 
fixed to 1 (a fully reflective surface). The overall normalization of 
the model is parameterized in terms of the luminosity of the seed 
black-body emitting surface at the thermal mound. In addition, by 
energy conservation the observed luminosity is equal to the grav¬ 
itational energy released in the accretion flow and proportional to 
the mass accretion rate. 

In order to account for emission from the polar cap outside 
the accretion flow, in addition to the F12 accretion flow spectral 
model (compmag) we include in our model one more black-body 
component (bb). The temperature of the black-body component 
describing the emission from the polar cap (kT^^) was initially 
tied to the temperature of the black-body seed photons from the 
thermal mound in the compmag model (kT™). However, we also 
tested if allowing the two black-body components to have different 
temperatures significantly improves the fit. This was the case for 
the spectra of 4U 0142-1-61 and RXJS 1708-4009, and the two pa¬ 
rameters were fitted independently. In the other 3 sources the two 
black-body temperatures were tied together. Finally, we included 
photoelectric absorption (phabs model in XSPEC) in order to ac¬ 
count for absorption by interstellar material along the line of sight. 
Therefore, our adopted model consists of two additive components 
seen through the same absorber: phabs * (bb + compmag). 

3.2 Spectral analysis 

The spectral analysis was performed with the XSPEC vl2.8 
spectral-fitting package (Arnaud et al. 1996), using the implemen¬ 
tation of the E12 model (model compmag). The spectral data were 
binned to have at least 20 net counts in each bin in order to allow 
the use of statistics. Before analyzing the spectra, we fitted them 
with the best-fit model published by Enoto et al. (2010) in order 
to check the consistency of our processing of the Suzaku data. In 
all cases we found consistent broad-band fluxes and spectral slopes 
above 10 keV with those reported in Enoto et al. (2010). 

In the case of IE 1841-045, which is located within a bright 
Supernova Remnant (SNR), we included in the spectral model a 
thermal plasma component (APEC; Smith et al. 2001) in order to 
account for the thermal emission from the SNR. In addition, the 
spectrum of 4U 0142-1-61 showed a strong artifact at an energy of 
1.84 keV (possibly due to the Si Ka absorption edge) that could 
not be modeled adequately. Eor this reason we excluded the 1.7- 
2.0 keV energy range from the spectral fit. 

The absorbed compmag + bb model gave good fits to the 
data (reduced < 1-2). The best fit parameters along with their 
90% confidence intervals for one interesting parameter (Ax^ = 
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2.709) are presented in Table|^ The unfolded spectra along with the 
best-fit models in E^{dN/dE) space are presented in Fig.[T] In the 
same figure we show the fit residuals in terms of the error of each 
data point (in a). These fits were performed with the parameter /3o 
(corresponding to the magnitude of the velocity of the accretion 
flow, in units of the speed of light, at the neutron star surface; Eq. 
13 fixed to /3o = 0 . 6 , as expected for free-falling material onto a neu¬ 
tron star of mass 1.4 Mq and radius 12.5 km. We also performed 
spectral fits with /3o free to vary. In all hut one case (IE 1547-5408) 
this parameter pegged at the maximum value (/3o ~ 1 . 0 ), while the 
hest-fit electron temperature kTe was unphysically low (~ 1 keV). 
This, in combination with the marginal improvement of the fit qual¬ 
ity over the fits with fixed fio, led us to adopt the best fit results with 
/3o fixed to 0.6, for all sources apart from IE 1547-5408 (however, 
fixing /3o to 0.6 in IE 1547-5408 also gives an acceptable fit). 

In the cases of 4U 0142-1-61 and IRXS 11708-4009 (see Fig.|3 
left), we see that the adopted model does not adequately reproduce 
the hard tail of the spectrum. In order to remedy this, we included 
in our model a thermal Comptonization component. Such a com¬ 
ponent could describe the Comptonization of the thermal emission 
from the neutron-star surface by an atmosphere of hot electrons 
(e.g. Ferrigno et al.. 2009). We used the model comptb in XSPEC 
(Farinelli et al. 2008; hereafter F08), with the parameter S fixed to 
0 (which corresponds to a pure thermal Comptonization spectrum). 
This model is parameterized in terms of the temperature of the seed 
hlack-hody photons which we allowed to differ from the 

temperature of seed black-body photons in the accretion flow, a 
power-law slope modifying the seed black-body spectrum (which 
we fixed to 7 = 3, i.e. a pure black body), the index of Green’s func¬ 
tion a, the temperature of the Comptonizing electrons (kTj^), and 
the relative contrihution of the seed black-body photons and upscat- 
tered photons to the total spectrum (the illuminating factor, log(A)). 
A spectral fit with a combination of the comptb, the compmag, 
and the bb models (model phabs* (comptb+compmag+bb) in 
XSPEC) greatly improved the fit, particularly at the hard part of the 
spectrum, while it gave more physical values to the other model pa¬ 
rameters. In order to limit the number of free parameters we fixed 
the terminal velocity of the accretion flow to 0 . 6 c, as with the spec¬ 
tral fits for the other sources. The model fit along with the residuals 
are shown in Eig. |3 (right), while the best fit parameters are pre¬ 
sented in Table |3 Hereafter we consider this as the best fit model 
for 4U 0142-1-61 and IRXS 11708-4009, and all derived parameters 
will be based on its best-fit parameters. 


4 DISCUSSION 

In the previous sections we have presented the results from the 
analysis of the X-ray spectra of 5 AXPs and SGRs with available 
data above lOkeV. The main driver of this investigation was to test 
whether the X-ray spectra of 5 AXPs/SGRs are consistent with the 
best available model for the X-ray spectrum expected from the ac¬ 
cretion column onto a pulsar (S ~ lO'^G; Earinelli et al. 2012). This 
model is better suited for the X-ray spectra of AXPs/SGRs than the 
models used in similar investigations in the past (e.g. Triimper et 
al. 2010; Guo et al. 2015), because of its more appropriate geom¬ 
etry and the use of (angle-averaged) magnetic cross-sections. Eur- 
thermore, it has a physical basis that enables us to derive and assess 
the validity of the physical parameters of the accretion flow within 
the context of the fail-back disk model, in contrast to phenomeno¬ 
logical fits using power-law models. 

Our analysis shows that the spectra of the 5 AXPs/SGRs stud¬ 


ied in this work are well fitted with a combination of a black-body 
model and the accretion column model of F12 (only in the cases 
of 4U 0142-1-61 and IRXS 11708-4009 we required an additional 
thermal Comptonization model). Even though this model includes 
several assumptions and simplifications in its treatment of the ac¬ 
cretion flow and radiative transfer (e.g. cylindrical geometry, angle 
averaged scattering magnetic cross sections, independent treatment 
of the thermal and bulk Comptonization) it is currently the best 
available model to study the X-ray emission from accreting pulsars 
and allows us to; (a) perform a basic test of the accretion model, and 
(b) obtain a first picture of the physical parameters of the accretion 
flow. 

In the cases of 4U 0142-1-61 and IRXS 11708-4009, the spec¬ 
tral fit requires an additional thermal Comptonization model, which 
we attribute to either the accretion flow below tbe shock or Comp¬ 
tonization of photons from the polar cap by electrons at its at¬ 
mosphere (c.f. Ferrigno et al. 2009). The low temperature of the 
seed photons (kT ~ 0.7 - 0.9 keV) and Comptonizing electrons 
(kTj^ ~ 0.1 -0.3keV) suggests that the Comptonization does not 
take place in the accretion flow (in contrast the electron tempera¬ 
ture in the accretion flow is ~ 22keV). However, we cannot rule 
out a contribution from the accretion flow. 

The analysis presented in §3 shows that this model can de¬ 
scribe very well the broad-band (0.5-200 keV) X-ray spectra of 
AXPs and SGRs, indicating that they could be produced from an 
accretion flow onto the pulsar. The spectral fits show that the tem¬ 
perature of the thermal emission from the polar cap is ~ 
0.2 - 0.6 keV. This is a factor of ~ 2 higher than the typical tem¬ 
peratures of “normal” pulsars, but it is consistent with the polar cap 
temperature of other AXPs (e.g. Aguilera et al. 2008), and could 
be understood in terms of heating of the polar cap by the fan beam 
(Triimper et al. 2013). 

4.1 Physical parameters of the accretion flow 

From the best-fit model parameters we can derive some basic in¬ 
formation about the physical conditions in the accretion column. 
Since the energy released is gravitational, we have that the mass 
accretion rate is given by 


where L is the luminosity of the emitted radiation, R and M are 
the radius and the mass of the neutron star respectively, and G is 
the gravitational constant. The mass accretion rates estimated from 
the luminosities listed in Table|4]and assuming a neutron-star mass 
M = 1.4M© and radius R = 12.5km are presented in Table|5] They 
are all in the range ~ lO'^g s“', which is ~ 100-1000 times lower 
than the accretion rates of typical accreting pulsars. They are also 
in good agreement with the fail-back disk model which predicts 
accretion rates of 10*^ - lO'^g s“* for ~ lO^yr after the supernova 
explosion (e.g. Ertan et al. 2009). 

The radius of the accretion column resulting form the spectral 
fits is generally between 160-330 m. Assuming that the accreted 
material is captured between the Alfven radius and the co-rotation 
radius, Kylafis et al. (2014) calculated that the radius of the accre¬ 
tion column is Vo ~ \lQM\f m, where Mu is the accretion 
rate in units of 10 *^ g s“*, and S 12 is the dipole magnetic field of the 
pulsar in units of lO'^ G. For the accretion rates given in Table |5] 
and a magnetic field of ~ lO'^ G, the above formula gives accre¬ 
tion column radii in excellent agreement with those estimated from 
the spectral fits. 
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Table 2. Spectral fit results 


Parameters 



Objects 




IE 1547-5408 

IE 1841-045+ 

4U 0142-H61 

IRXS 11708-4009 

SGR1900-tl4 

Nh (lO^^cm-^) 

9 qi +0.12 
^•^^-0.08 

2.66±0.1 

0.68 ±0.01 

1.61 ±0.06 

2.6 ±0.4 

krPC (keV) 

0.61 ±0.03 

0.59 ±0.01 

0.275 ±0.004 

0.197 ±0.005 

0 41+012 

^•‘ti-0,09 

Normbb (xIQ-^Lss/D^) 

4.76±0.50 

3.01 ±0.02 

15 9^-^ 

9.9±0.3 

u./O-0 29 

kT™ (keV) 

0.61* 

0.59* 


0.397 ±0.006 

0.41* 

kXA (keV) 

0.23!^;® 

60^1 

-0.001 

0.02±0.01 

0.2(< 31.7) 

Tj-t 

170+60 

i /^_10 

156!^* 

104+*’ 

160 ±3 

^^•^-60 

/3o (c) 

0 47-^<^-02 
^•^'-0.08 

0.6(f) 

0.6(f) 

0.6(f) 

0.6(f) 

ro (m) 

272^ 

28000 250) 

169 ±2 

183.4±0.4 

253+76 

^^•^-60 

Normcom pmag 

321+^8 

-192 

yo_23 

17600!’™ 

7227 ±51 

1 o 1 ±288 

X^ldof 

1398.7/1341 

1112.9/962 

1034/919 

1251/1078 

85.2/85 


t The model for the spectral fit of IE 1841-045 includes an APEC thermal plasma component (e.g. Smith et al. 2001; Foster et al. 2012). Its best-fit 
temperature and abundance are kT = 0.51 zb 0.02 keV, and Z = 0.40 di 0.07 respectively (c.f. FigD. 

* The temperature of the black-body component is fixed to the temperature of the black-body parameter of the compmag component. 

% The temperature of the electrons in the accretion flow; see §4.2 for a detailed discussion of the fit results. 

I The Thompson optical depth along the accretion flow (parallel to the magnetic field) (with a conversion factor of (j| | = 10“^err). 

** The normalization of the compmag model is defined as where is the radius of the black-body emitting area, and D is the distance of the 

source in units of lOkpe. 

(f) Fixed pai'ameter. 


Table 3. Spectral fit results for 4U 0142+61 and IRXS J1708-4009 


Parameter 

Eit 

results 


4U0142-t61 

IRXS 11708-4009 

Nh (lO’+m-’) 

0.67 ±0.01 

-0.13 

kT-TC (keV) 

0.28^:<>t 

^•^^+0.02 

NomiTc (x 10'^T39/D^) 

29.6!«;’ 

22 7'^-^ 
^^•'-0.1 

TCtt 

'bb 

3.0 

3.0 

TC 

9 99+0.14 

^•^'-0.69 

1 4-^0.3 
^•^-0.2 


0.0 

0.0 

log(A)’''* 


0 

b 

c^ 

krjet (keV) 

0.88^:>^ 

0 74‘^'i* 

-0.08 

q-jt 

co±356 

jj_23 

30+35 

-27 

/3o (c)^^ 

0.6 

0.6 

kT™ (keV) 

^ -^^-0.09 


kTA (keV) 

22 ;’ 

40.2+™ 

fo (m) 

329 ±9 

231!^*’ 

N Ornicom pmag 

1 7+18.8 
^•'-0.13 

7 9+22.6 
'•^-4.4 

X^/dof 

1001/919 

1258/1106 


11 Model parameter fixed. 

f The temperature of the electrons in the comptb model. 

J The Thompson optical depth along the accretion flow (parallel to the magnetic field) (with a conversion factor of (T|| = Kr^erj). 

1 The temperature of the electrons in the accretion flow; see §4.2 for a detailed discussion of the fit results. 

* The normalization of the compmag model is defined as where is the radius of the black-body emitting area, and D is the distance of the 

source in units of 10 kpc. 

Key for the formation of the part of the spectrum extending n (z)- ^ (4) 

above lOkeV is, as discussed in §3.1, a transverse optical depth of v(z)nrgmp 

at least ~ 1. From the continuity equation we have that the electron 
density ne(z) at a given distance z from the center of the neutron 

star is (c f Eq 43 of FI2)' where M is the mass infall rate, v(z) is the velocity of the accretion 

flow at distance from the neutron star center z, n is the radius of 
the accretion column, and nip is the proton mass. 
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Figure 1. The Suzaku (XIS and HXD) and INTEGRAL (ISGRI) spectra of IE 1547-5408, IE 1841-045, and SGR1900+14, (in E^^ scale) along with their 
best-fit models presented in Table[2| The Suzaku XIS and HXD spectra are shown in black and red respectively, while the INTEGRAL ISGRI spectra are shown 
in green. For clarity we only show spectra from one XIS detector. In the case of IE 1547-5408 the quality of the INTEGRAL ISGRI spectra was poor and only 
the Suzaku HXD-PIN spectra ai‘e shown. The contribution of the black-body component representing the polar cap emission, and the compmag accretion flow 
component to the overall spectral model (solid line), are shown by the dashed lines (the spectnim of IE 1841-045 also includes a thermal plasma component 
shown by the dash-dot line). The bottom panel of each figure shows the residuals from the best-fit model (in a). 


Table 4. Source Luminosity 


Object 

Distance^ 


Flux 

Luminosity 



(0.1-10.0) keV 

(10.0-250.0) keV 

(0.1-10.0) keV 

(10.0-250.0) keV 



obs. (unabs.) 

obs.^ 

obs. (unabs.) 

obs. 


kpc 

(10“^* erg cm“^ s~^ 

') (10“*'erg ) 

(10**^ erg s“*) 

(10**^ erg s“*) 

IE 1547-5408 

4.5 

6.67 (12.66) 

19.2 

1.6 (3.1) 

4.6 

IE 1841-04* 

8.5 

2.13 (4.43) 

8.52 

1.8 (3.8) 

11.0 

4U0142+61 

3.6 

13.8 (31.7) 

1.3 

2.1 (4.9) 

0.2 

IRXS J1708-4009 

3.8 

4.12(26.3) 

8.2 

0.7 (4.5) 

1.4 

SGR1900+14 

12.5 

0.54(1.4) 

1.86 

1.0 (2.6) 

3.5 


t The adopted distances are from Olausen & Kaspi, 2014. 
f The unabsorbed hard X-ray flux is equal to the absorbed flux for the Hi column densities considered here. 

* The reported fluxes for IE 1841-04 exclude the thermal emission from the suiTounding Supernova Remnant. 


Using the velocity profile given in Eq. [3 the above equation 
becomes 


ne{z) = 


M 

TvcmprlPo 


lyr, 


(5) 


center of the neutron star R), and rj is the exponent of the velocity 
profile (in our case fixed to 0.5). 

Then the transverse optical depth at distance z is 


T± = Heroar = 


LR 

'KCwipGMraPa 



( 6 ) 


where c is the speed of light, /3o is the velocity (in terms of the 
speed of light) at the surface of the neutron star (distance from the 


Based on the hest-fit parameters given in Table we find 
that the transverse optical depths at the base of the accretion col- 
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Figure 2. The Suzaku (XIS and HXD) and INTEGRAL (ISGRI) spectra of 4U0142+61 and IRXS J1708-4009 (in scale) fitted with the phabs* (bb + 
compmag) (left; Tabled and the phabs * (comptb + compmag) (right; Table[3j models. The Suzaku XIS and HXD spectra are shown in black and red 
respectively, while the INTEGRAL ISGRI spectra are shown in green. For clarity we only show spectra from one XIS detector. The spectra have been corrected 
for instrumental effects. We clearly see how the addition of the thermal Comptonization model improves the fit. The bottom panel of each figure shows the 
residuals from the best-fit model (in a). 


Table 5. Physical parameters of the accretion flow 


Object 

Dist.f 

kpc 

M 

lO'^grs^' 

'TJ. 

Rpc 

km 

IE 1547-5408 

4.5 

5.2 

1.4 

0.77 

IE 1841-04 

8.5 

10.0 

0.25 

3.6 

4U 0142-1-61 

3.6 

3.4 

0.7 

18.3^ 

IRXS J1708-4009 

3.8 

3.9 

1.2 

62.3 

SGR1900-I-14 

12.5 

4.4 

1.2 

5.4 


t The adopted distances are from Olausen & Kaspi, 2014. 

^ The radius of the polar cap for 4U 0142-1-61 and IRXS J1708-4009 is based on the best-fit normalization and black-body temperature of the comptb 
model and it is likely overestimated since this model may also account for some emission produced in the accretion flow. 


umn are 1 (Tabled. Kylafis et al. (2014) demonstrated that 
even in such relatively low optical depths, thermal and bulk-motion 
Comptonization can produce a hard spectrum with a photon index 
r ~ 1-1.5 extending up to ~ 200 keV. On the other hand, the best- 
fit Thompson vertical optical depth along the accretion flow is typ¬ 
ically r ~ 100-400. 


4.2 Limitations of the model 

Although the accretion column model provides good fits to the 
spectra of the 5 AXPs/SGRs studied in this work, some of the 


model parameters take unphysical values. The two most extreme 
examples are the magnitude of the velocity, in units of the veloc¬ 
ity of light, that the gas in the accretion flow would aquire at the 
neutron star surface (/3o), which reaches unity if left as a free pa¬ 
rameter, and the electron temperature in the accretion flow {kTZ), 
which is much lower than 1 keV (Tables |2[^. This indicates that 
the spectral model produces the high-energy photons only through 
the bulk-motion Comptonization mechanism, while the contribu¬ 
tion of thermal Comptonization is negligible, even when the energy 
available through the bulk-motion Comptonization mechanism is 
not sufficient to produce the hard X-ray photons (as in the case 
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when /3o was fixed to a lower value; see §3). This behavior could 
be the result of the following factors: 

(a) the seed photon spectrum assumed in the F12 model is a 
blackbody produced at the base of the accretion column, and it ig¬ 
nores the contribution of the bremsstrahlung radiation produced by 
the electrons below the shock. In fact, the F12 model does not take 
into account at all the flow below the radiative shock, but instead 
assumes that the accretion flow (above the shock) has a uniform 
temperature Tg. For consistency, the model ought to take into ac¬ 
count the bremsstrahlung radiation that is produced from the elec¬ 
trons in the accretion flow (above the shock), where the same elec¬ 
trons upscatter the produced photons, thus self-consistently linking 
the observed spectrum below ~ 20 keV with the hard X-ray spec¬ 
trum. Since a bremsstrahlung component would contribute in terms 
of a harder seed spectrum, the required contribution of bulk-motion 
Comptonization would be reduced resulting in more realistic values 
for the terminal velocity of the accretion flow. 

(b) the electron temperature is decoupled from the parame¬ 
ters of the accretion flow. In a more realistic scenario one would 
expect that the same electrons that are participating in the bulk- 
motion above the shock are also responsible for the thermal Comp¬ 
tonization below the shock (Lyubarskii & Sunyaev 1982; Kylafis 
et al. 2014), resulting in much higher electron temperatures and 
harder Comptonized spectra, reducing the need to interpret the en¬ 
tire hard spectrum through bulk-motion Comptonization. 

(c) the mass accretion rate estimated through Eq. (3) as¬ 
sumes isotropic emission, which obviously would be inappropriate 
if beaming effects are coming into play. In their study of the en¬ 
ergy dependent pulse-profiles of 4U 0142-1-61, the brightest AXP, 
Triimper et al. (2013) estimate a luminosity for the fan beam of 
1.8 X 10^^ erg s“*, and for the polar beam 0.3 x 10^^ erg s“', re¬ 
sulting in a total luminosity of 2.1 x 10^^ erg s^'in the 0.8-160 keV 
energy range. This is 40% lower than the luminosity estimated 
from Eq. (3), indicating that the neglect of beaming effects may 
significantly overestimate the accretion rate, and subsequently any 
parameters depending on them (e.g. optical depth). 

Nonetheless, this model is the best available spectral model 
for the accretion flow in an accreting pulsar, and it provides a good 
representation of the observed spectra. Improved models address¬ 
ing the above limitations would help to obtain more precise infor¬ 
mation for the parameters of the accretion flow, but they should not 
change the general picture presented here. 

4.3 AXPs and SGRs as accretion-powered pulsars 

The debate regarding the interpretation of the quiescent emission of 
AXPs and SGRs as magnetars would largely benefit from the mea¬ 
surement of their dipole megnetic field strengths from cyclotron 
lines. For example, more sensitive observations of IE 1841-04, 
which show evidence for an absorption feature at the 3a level (An 
et al. 2013), can give insights into their nature by comparing their 
spectral and timing properties with those expected from different 
models. In addition, recently, Tiengo et al. (2013) reported an in¬ 
ferred magnetic field of > 2 x lO'"* G in SGR 0418-1-5729 assuming 
a proton cyclotron line, which is interpreted as a multipole mag¬ 
netic field; on the other hand, if an electron cyclotron line is as¬ 
sumed, the strength of the magnetic field is closer to ~ lO" G. 
The fact that the broad-band X-ray spectra of five AXPs/SGRs with 
hard X-ray emission can be well reproduced by an accretion flow 
model onto a pulsar, and the resulting parameters are consistent 
with those expected from the fail-back disk model, despite the lim¬ 
itations of the spectral models, gives further support to the fail-back 


disk model. Eurthermore, fits with the original version of the BW07 
model, which includes bremsstrahlung seed photons (but has strong 
correlations and degeneracies between the model parameters) also 
provide an equally good description to the spectral data. 

Triimper et al. (2013) find that the energy-dependent pulse 
profiles of the brightest AXP, 4U 0142-1-61, can be explained by 
the emission of a fan beam originating in the accretion shock and 
a polar beam from the polar cap. The latter is produced by the il¬ 
lumination of the polar cap by the fan beam, which is subject to 
gravitational bending towards the neutron-star surface. This model 
allows us to determine the inclination of the spin axis with respect 
to the line-of-sight (i = 60°), the angle between the magnetic and 
spin axes (a = 30°), and the height of the accretion shock from 
the neutron star surface {H = 2 km, corresponding to z = 14.5 km 
from the neutron star center). The combination of the constraints 
from the timing analysis, with the results from the energy spectra 
presented in this work, provide support to the accretion model for 
AXPs/SGRs from two independent, but complementary perspec¬ 
tives. Development of more sophisticated spectral models for the 
accretion flow as well as modeling of the phase resolved spectra 
(e.g. Triimper et al. 2013) for all sources with good quality hard 
X-ray spectra will allow us to further test the fall-back disk model 
for this class of objects and better constrain the parameters of the 
accretion flow. 

5 CONCLUSIONS 

Our main conclusions are enumerated below. 

(i) We presented a systematic analysis of the X-ray spectra of 5 
AXPs/SGRs. We fitted their Suzaku and INTEGRAL data with the 
spectral model of P12 for the radiation from material infaling onto 
a pulsar through an accretion column. We found that this model 
gives good fits to the observed broad band (0.5-200 keV) spectra, 
supporting the fail-back disk model. 

(ii) We measured accretion rates in the ~ lO'^gs^* range, also 
in agreement with the calculations of Ertan et al. (2009). 

(iii) We found that the Thompson optical depth along the accre¬ 
tion column is typically ^ 100-400, while the transverse optical 
depth is ~ 1.0. 

(iv) The radius of the accretion column is ^ 200 m, consistent 
with the estimation of Kylafis et al. (2014) for the fall-back disk 
model. 

(v) These results are in very good agreement with the indepen¬ 
dent constraints on the X-ray beaming accretion flow geometry 
from the analysis of the energy-dependent pulse profiles (Triimper 
et al. 2013). 

(vi) In general, these results show that the observed broad-band 
X-ray spectra of AXPs and SGRs showing significant X-ray emis¬ 
sion above lOkeV are consistent with the fail-back disk model. 
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